Long-time properties of MHD turbulence and the role of symmetries 
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We investigate long-time properties of three-dimensional MHD turbulence in the absence of forcing 
and examine in particular the role played by the quadratic invariants of the system and by the 
symmetries of the initial configurations. We observe that, when sufficient accuracy is used, initial 
conditions with a high degree of symmetries, as in the absence of helicity, do not travel through 
parameter space over time whereas by perturbing these solutions either explicitly or implicitly using 
for example single precision for long times, the flows depart from their original behavior and can 
become either strongly helical, or have a strong alignment between the velocity and the magnetic 
field. When the symmetries are broken, the flows evolve towards different end states, as predicted 
by statistical arguments for non-dissipative systems with the addition of an energy minimization 
principle, as already analyzed in J] for random initial conditions using a moderate number of Fourier 
modes. Furthermore, the alignment properties of these flows, between velocity, vorticity, magnetic 
potential, induction and current, correspond to the dominance of two main regimes, one helically 
dominated and one in quasi-equipartition of kinetic and magnetic energy. We also contrast the 
scaling of the ratio of magnetic energy to kinetic energy as a function of wavenumber to the ratio 
of eddy turn-over time to Alfven time as a function of wavenumber. We find that the former ratio 
is constant with an approximate equipartition for scales smaller than the largest scale of the flow 
whereas the ratio of time scales increases with increasing wavenumber. 

PACS numbers: 47.65.-d , 47.27.Gs 47.27.ek , 94.05.Lk 

I. INTRODUCTION 
A. The context 

Magnetic fields pervade the universe and often play an essential role in configuring and constraining structures, as in 
the case of intergalactic jets, or closer to Earth, in the Solar Wind or the magnetosphere. Magnetic pressure contributes 
to the containment of the heliosphere [2], and it may retard, together with turbulent pressure, the gravitational collapse 
of molecular clouds in the interstellar medium. Magnetic fields are also known to accelerate the motion of charged 
particles in the magnetospheres of the planets in the Solar System, as for example in the case of Jupiter's aurora [5] 
(see [1] for the main characteristic of the auroral emissions as of today). 

Such magnetic fields have been observed in a variety of media to be turbulent, such as in the Solar Wind [5] (see 
[B] for a recent review) , in the interstellar medium where it is thought to be responsible for strong velocity shear and 
intermittency [7J, or more recently in the heliosheath [§]. Using CLUSTER data with short separation between the 
satellites, it was shown in [9] that anisotropy of the energy Fourier spectra develops at small scales as predicted in 
weak MHD turbulence theory [lOHE]. Solar Wind turbulence can also help focus Langmuir wave packets which are 
routinely observed using Ulysses or STEREO spacecrafts |13j . 

Furthermore, magnetic fields can lead to extreme energetic events due to reconnection of magnetic field lines in 
highly turbulent media; solar flares are one such example, the prediction of which is one of the purpose of space 
weather research because of the disturbance to Earth's communication networks and power grids. The penetration 
of the Solar Wind into the Earth's magnetosphere can be explained by the development of Kelvin-Helmoltz (KH) 
vortices, as observed |14j . and a relationship between such KH instabilities and flux transfer events was found recently 
during substorms using multiple spacecrafts |15j . 

An understanding of both fluid and MHD turbulence has escaped us for a long time. Is MHD turbulence similar 
to hydrodynamic turbulence, with a Kolmogorov energy spectrum, Exuik) ~ e 2 / 3 fc -5 / 3 (hereafter K41), with e = 
dE / DT the energy dissipation rate, perhaps with an anisotropy due to the presence of strong uniform fields of 
magnitude Bq? Or is it different, because Alfven waves propagate that alter and dampen the nonlinear dynamics 
of turbulent flows, leading to a so-called Irsohnnikov-Kraichnan energy spectrum, Ejx{k) ~ [eBo\ 1 / 2 k~ 3 '/ 2 (hereafter 
IK)? And is there one answer to these questions, or is universality broken in MHD, as sometimes advocated? For 
example, it was found in [TB] that one can observe three different energy spectra (K41, IK and weak turbulence, 
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hereafter WT, Ewr(k±) ~ fcj 2 ), for three different initial conditions of the magnetic field, using the same velocity 
field, and with the same ideal invariants, namely total energy, total magnetic helicity and total cross correlation 
between the velocity and the magnetic field, with moreover no imposed external field, no forcing, unit magnetic 
Prandtl number, and equal kinetic and magnetic energy initially? In other words, for such a set of initial conditions, 
nothing allows for distinguishing these three configurations from a statistical point of view, no externally imposed 
time-scale is present, and the only constraint is that numerically the flows and fields are forced to follow the four-fold 
symmetry of the initial conditions. Using this symmetry, one can gain in resolution and cost of computation, and 
thus the Reynolds numbers are quite large (Taylor Reynolds number in excess of 1200), with equivalent resolutions of 
2048 3 grid points. Similar results are found to hold in the forced case as well, for which long-time averaging is feasible 

Em. 

The difference between these three power laws could be due to non-local interactions in Fourier space, between 
widely separated scales. Non-local interactions are thought to be more prevalent in magnetohydrodynamics (MHD) 
than in hydrodynamics, as measured in high-resolution numerical simulations |18Lll9j . and such non-locality in Fourier 
space is advocated in the differentiation between a K41 and an IK spectrum, but what would make one of the three 
flows studied in |16) more non-local than others? Perhaps the different behaviors come from another factor. On the 
one hand, it could be that the ratio of kinetic to magnetic energy, in particular in the gravest mode, matters, as 
indicated in [16]. On the other hand, the invariants, which are quadratic in the basic fields, are identical but higher 
order moments could differ; for example, it was shown in [20] that the skewness of one of the flows studied in [16] is 
measurably larger than that of the two other flows, when looking at both the velocity and the magnetic field. 

The assumption that with the same invariants, the three flows should behave in similar ways is based on an 
assumption of ergodicity. However, the ergodicity of turbulent flows has been put into question in a variety of contexts. 
It has been observed that long-time memory effects can be found in such flows, for example in two-dimensional (2D) 
MHD turbulence |21j , where large bursts of energy were observed to evolve on time scales of the order of one hundred 
turn-over times; it was also found more recently in numerical simulations of 3D hydrodynamic (HD) turbulence |22] 
and in laboratory experiments (see e.g. [23] and references therein). In the atmospheric boundary layer, one observes 
that statistics can depend on the large scales; this may be related to averaging over regions with local fluctuations 
in Reynolds number [23]. Such a transfer between non-local triads is shown to lead, however, to local exchanges of 
energy in hydrodynamics |25) . 

The persistence of modes for long times in MHD, associated with large-scale coherent structures, was recently linked 
to normal modes appearing because of magnetic helicity, a large-scale invariant in the ideal case, and leading to an 
apparent breaking of ergodicity at the largest scale, insofar as these structures persist for long times |26j . A similar 
phenomenon occurs for hydrodynamics in the presence of solid body rotation: although the Coriolis force is linear, 
it affects the dynamics of rotating turbulence in slowing it down substantially |27j ; this can be attributed to inertial 
waves, nonlinear transfer occurring only through (quasi-) resonances |12j . 

In selective decay, some invariants are viewed as more sturdy than others; thus, they may influence the long-term 
dynamics of decaying turbulent flows. This hypothesis is based on the fact that invariants may have different physical 
dimensions, for example magnetic potential < A 2 >, with b = V x A, and total energy | < |v| 2 + |b| 2 > in 2D, or 
magnetic helicity < a • b > and total energy in 3D: since dissipation involves a Laplacian, it is thought that < A 2 > 
or < Hm > will decay more slowly than energy. However, the third invariant in ideal MHD, He =< v • b >, has the 
same dimension as energy and it could also influence the long term dynamics, becoming strong in relative terms, that 
is with respect to the energy, implying an alignment between the velocity and the magnetic field, a phenomenon called 
dynamic alignment. The relative importance of these two effects was explored, both theoretically and numerically in 
[T] [25] . These theoretical considerations based on statistical mechanics of a truncated system of modes were backed up 
by rather low resolution numerical simulations which nevertheless clearly demonstrated the validity of the approach: 
the end state of such flows was determined by the respective ratio of their three invariants. Will the same happen 
here, when starting with the three flows studied in [16] . which statistically are equivalent but which display different 
inertial range dynamics at peak of dissipation (and in the statistically steady state as well)? This is the main question 
that this paper is addressing, using direct numerical simulations of the MHD equations in three space dimensions. 



B. The equations 



We now give the MHD equations for an incompressible fluid with v and b respectively the velocity and magnetic 
fields in Alfvenic units: 

dv 1 

v+v-Vv= W + j x b + zA7 2 v, (1) 

at po 

f)h 

— = V x (v x b) + r;V 2 b; (2) 
at 
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po = 1 is the (uniform) density (and b is then dimensionally a velocity as well, the Alfven velocity), V is the total 
pressure, V ■ v = V ■ b = 0, and v and r\ are respectively the kinematic viscosity and magnetic diffusivity; we take 
v = ij. With v = 0, r) — 0, the energy Et, the cross helicity Hq and the magnetic helicity Hm, defined as 

E T = E V + E m = (v 2 + b 2 ) /2 , H c = (w b) /2 , H M = (A ■ b) /2 , 

are conserved. Relative helicities can be defined as follows: 

p c = cos [v, b] , p M = cos [A, b] , p v = cos [v, w] ; (3) 

they correspond to the degree of alignment between various vectors: the velocity, the magnetic field, the magnetic 
potential (with pu = ±1 defining a force-free field), or the vorticity (with pv = ±1 defining the so-called Beltrami 
configuration). In the latter case, the relative kinetic helicity involves the vorticity, uj = V x v; the total kinetic 
helicity is an invariant of the Euler equations (b = 0, v = 0). 

The kinetic energy spectrum is the Fourier transform of the velocity two-point correlation function. Once homo- 
geneity, isotropy and incompressibility have been taken into account, only two defining functions remain: Ey(k) is 
proportional to the kinetic energy, with J Ey(k)dk — Ey = ^(v 2 ), and the kinetic helicity, Hy(k), stems from the 
anti-symmetric part of the velocity gradient tensor. Similar definitions hold for the magnetic and cross correlation 
functions (note that helicity is a pseudo-scalar). Finally, the kinetic and magnetic Reynolds numbers are defined as 

Ry = U L Q /v . R M = U L /r), 

where Uo, Lo are the characteristic velocity and length scale. The integral scale is defined as 

_ J[Ey(k)/k] dk 
" lt ~ jEy(k)dk ■ 



C. Predictions from statistical mechanics 



The statistical equilibria in 3D MHD were derived in [21] . They are the long-time solutions to a truncated system 
of Fourier modes, with k m i n and k max the minimum and maximum wave numbers respectively; these modes are 
coupled through the nonlinear ideal MHD equations (v =0, r] = 0), and subject to the conservation of all quadratic 
invariants. Defining a / 0, /5 and 7 as the Lagrange multipliers associated with the Et, Hm and He invariants, 
namely uEt + (3Hm + 7-ffc , these equilibria read, assuming that the magnetic helicity is non-zero (f3 ^ 0) : 

Hm ^ = -^ Wk) 5 H j( k ) = k2H M( k ) ; H c (k) = ^- Hj(k) , H v (k) = ^ Hj(k) , (4) 



Eu(k) = H j{ k) = , Ey(k) = (r 2 V(k) + E M (k) = (l - E M (k) , (5) 

where Hj — J k 2 HM(k)dk is the current helicity, and 

«>o, r 2 = i-^>o ; v{k)= (i-^ip) >° > vfe e [k mi n,k max } ■ (6) 

Hm and He are not definite positive, and furthermore, Hm does not have the same physical dimension as Et and 
He, and hence ft does not have the same physical dimension as a and 7. In order to fulfill realizability conditions 
(positivity of energy, and Schwarz inequalities involving the helicities), necessary relationships between coefficients 
can be derived, involving k m i n (see [2§]). 

When (3 = and thus ffjvf(fe) = 0, Hj(k) = 0, one finds that the kinetic helicity is also equal to zero, that we have 
equipartition of energy at all wave numbers with EM(k) = Ey(k) = 8irk 2 /(aF 2 ), and that Hc(k) = —4:TT^k 2 /(a 2 r 2 ); 
thus the relative cross helicity 2Hc(k) / Ey{k) is constant in that case. When 7 = and thus Hc(k) = 0, the 
kinetic helicity is also equal to zero and the kinetic energy has its non-helical expression, Ey{k) = 8irk 2 /a; the 
magnetic energy and helicity can peak at low wavenumber when j3 is large enough, and the relative magnetic helicity 
kHM(k)/ Ej^ik) ~ 1/fc, i.e. it is stronger in the largest scales of the flow, a result that persists in the general case 
08^ 0,7^0). 
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When considering Hj(k) (instead of Hj^(k)), note that all Fourier spectra are strictly proportional, with coefficients 
uniquely determined by initial conditions given the values of the invariants, except for the kinetic energy; also note 
that, Vfc, one has Ey(k) < iJjyf(fc), the equality arising only when there is either no magnetic helicity, or maximal 
cross-correlation, or for k max — > oo. Similarly, the residual helicity defined as 

H R (k) = H v (k) - Hj(k) = -T 2 Hj(k) , 

is of the sign opposite to that of the current and of the magnetic helicity, and Hn(k) becomes equal to zero only 
for maximal cross-correlation (r 2 = 0), except for the trivial non-helical case of course. The relative helicity H^(k), 
integrated over the small scales, is the motor of the nonlinear dynamo problem, i.e. the growth of large-scale magnetic 
energy because of small-scale helical motions; note that Hr reduces to the kinetic helicity in the kinematic regime 
when the magnetic field is weak, thus recovering the so-called "alpha" effect (see [30] for a recent comprehensive 
review) . 

It was shown in [TJ [2H] that these solutions can be seen as indicators of the long-time behavior of 3D MHD systems 
left to decay, because of a principle of minimization of total energy. Three main regions of parameter space can be 
seen as attractors to the dynamics: a magnetic helicity dominated region, an alignment (strong He) region, and an 
intermediate region. The relaxation principle is well founded when there is magnetic helicity in the system, since 
dimcnsionally Hm weighs the large scales more so than the energy or the cross-correlation, but when Hm = 0, it is 
not so clear what happens. It can be shown, using minimum energy principles following [31 , that the resulting fields 
are u = and j cx b when there is no cross-correlation (the constraint is simply that Hm remain constant) , whereas in 
the general case, the solution is a bit more involved (see eq. (3.13) sq. in [35]; see also 28 J. The main purpose of this 
work is to investigate this long-term dynamics when considering the three initial conditions used in [16) which have 
the same quadratic invariants (Et = 1/4, Hm = and Hq < 4% in relative terms) and thus presumably the same 
final asymptotic state and yet, at peak of dissipation, show clear differences in their inertial range scaling. We shall 
also investigate other deterministic flows with either cross-helicity or magnetic helicity to see whether they evolve as 
well towards these attractors. 



D. Description of the initial conditions for all the computations 

Table [I] summarizes the main characteristics of the 36 runs described in this paper. Further details on the com- 
putations are given below, when specifying the initial velocity and magnetic fields. All runs use the Geophysical 
High-Order Suite for Turbulence code (GHOST) [33J unless otherwise stated in the "Remarks" column; TYGRS 
stands for a code which implements the four-fold symmetries of the Taylor-Green (TG) flow and its extensions to 
MHD [16]. Values at t = T wo = 100t nl of 2H C /E T , H M /E T , and E T , with t nl = L /U , are given in the Table; 
note that k m i n = 1 has been taken as a normalizing factor for the ratio involving magnetic helicity. For all runs v = t] 
and Em = Ey = 0.125 initially. The groups divided by horizontal lines correspond to different line styles in Figs. [T] 
[2J and|4| All runs are performed on a grid of 64 3 points with 1/ = 2 x 10~ 3 , except for runs RIOc and Rllb done on 
grids of 32 3 points, runs RIOd, R17b, and R19 done on grids of 128 3 points, run R21 done on a grid of 192 3 points, 
and runs R9b and R23b done on grids of 256 3 points (see the last column for the value of the viscosity in these cases) . 
/D and /S stand for double and single precision. I, A and C in the second column refer to the Taylor-Green flows 
studied in [TB]: I, A and C are for the insulating boundary conditions (I), the alternate insulating conditions (A) 
and the conducting one (C); by insulating or conducting it is meant that in the box in which the computations are 
performed, the current is either parallel or normal to the walls. In most cases (except those labeled "TYGRS"), the 
symmetries of the Taylor-Green initial conditions are not enforced and can be broken. 

Runs in which noise of amplitude 10~ x relative to the host flow has been added to both the kinetic and magnetic 
energy are denoted "+10 _a " in the second column. Furthermore, "+V=xABC"/ or "+B=xABC"/ or "V&B=xABC" 
indicate that a Beltrami ABC flow has been added to either the velocity, magnetic field, or both for the initial 
conditions, with "x" indicating the fraction (in terms of energy) of the initial condition which is ABC. "TGx;ABCy" 
stands for a modified Taylor-Green velocity at k — x and a magnetic field which is a Beltrami ABC flow at k = y. OT 
stands for the Orszag-Tang vortex generalized to three dimensions as the initial condition studied in [34], and finally 
the last 4 runs have the velocity of [33] and a mixture of OT and ABC with the specified fractions for the magnetic 
field. The purpose here is to be able to vary the cross-helicity He and the magnetic helicity Hm of well-studied 
configurations in MHD turbulence, at a fixed total energy Et, the same in all runs; indeed, all computations have 
equal initial kinetic and magnetic energy, with Et = Ey + Em = 1/4. 

GHOST is a general purpose pseudo-spectral community code with periodic boundary conditions; the code is now 
parallelized up to ~ 98, 000 processors, using a hybrid (MPI-Open-MP) methodology that becomes advantageous at 
high resolution [33]. Runs Rib, R5b, and RlOb are done using a similar code, TYGRS, but in which the four-fold 
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symmetries of the Taylor-Green configuration are enforced at all times [35., ,36]; TYGRS follows the same parallelization 
methodology as GHOST. The Taylor-Green velocity is: 

v x — v sin/c^x cosk v y cos k v z 7 v y — — v cosk v x smk v y cosk v z, v z =0 , 

and the three different initial conditions for the magnetic field are in that case: 

b\. = bQCOsk m x sink m y sinfc m z, by — &gSinfc m :E cos k m y sinfc m z, b\ = —2b\ ) s\~n.k m x smk m y cos k m z ; 

l>x = b^cosk m x sinfc m y sinfc m z, by — — b A sink m x cos k m y smk m z, b A = ; 

and 

b G = b c sin k m x cos k m y cos k m z, by = b° cos k m X smk m y cos k rn z, b c = — 2b c cos k m x cos k m y sinfc m z . 

When computations in which these initial fields are perturbed with an added noise, the amplitude of that noise 
relative to the energy in the Taylor-Green initial condition is indicated in the second column. This noise has randomly 
generated phases with an energy spectrum of the form: 

AT AT ( (logfc-logfco) 2 

N = N exp — ^ 

V 2 (logo-) 2 

In all cases the noise is centered around fc = 2 and has a = 2. Noise of this form is added to both the magnetic 
and kinetic energy and introduces small perturbations in the initial magnetic helicity and cross helicity relative to the 
total energy depending on the random phases generated and the amplitude of the noise. 

We also performed some runs which have significant amounts of helicity since helicity is a main indicator of 
the behavior of such flows, at least in the ideal regime. Since the Taylor-Green runs have no helicity, different 
configurations are also studied. The "TGx;ABCy" type is one for which the velocity is a modified Taylor-Green, such 
that wJ G = —Vy G , Vy G — — «J G , and v^ G — 0, centered at wavenumber k v — x with x equal to either 2 or 3, and 
the magnetic field is a Beltrami ABC flow centered at wavenumber k m — y with y equal to either 1, 2, or 3; the ABC 
magnetic field is: 



K L = &o L [5cos(fc m y) + Csin(fc m z)] . 
K BC = b o BC [Asm(k m x) + C cos(k m z)} , 
bf BC = ^ sc [Acos(fc m x) + £sin(/c m y)] . 

The OT configuration is that of the generalization of the Orszag-Tang vortex to three dimensions, as studied in 
[31], with the velocity and magnetic fields defined as: 

v° T = -2v^ T sin k v y, v° T = 2v JT sin k v x, v OT = , 

and 

b GT = b GT [—2a sin 2k m y + a sin k m z] , b y T = b T [2a sin k m x + a sin k m z] , b GT = b T [a sin k m x + a sin k m y] ; 

the parameter a allows one to modify the cross-correlation between the two fields; the choice a = 0.8 gives a relative 
correlation of 0.41. 

Finally, initial conditions which are mixtures of the above types are also studied. Runs R15 and R16 have the 
Taylor-Green velocity and perturb the Taylor-Green magnetic field initial condition with an "A" configuration, with 
an ABC Beltrami field such that: 

b A+ABC = IjA+abc (£ lC0S k ml x sinfc m iy sin k„ a z + £, 2 [B cos(k m2 y) + C sm(k m2 z)}) , 
frA+ABC = b A+ABC (-^smkmxx cosk ml y smk m iz + £ 2 [Asm(k m 2x) + Ccos(fc TO2 2)]) , 



b A+ABC = b A+ABC^ [AcQs{km2x) + Bfjin{km2y)] 



The parameters £i and £2 set the relative fractions of the Taylor-Green and ABC portions of the initial condition. 
Run Rf5 has £1 = 0.99 and £2 = 0.01 and run R16 has £1 = 0.9 and £2 = 0.1. Both of these flows have k m \ and k m2 
such that both portions of the initial condition are at k = 3. Initial conditions such as these allow for a perturbation 
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in the magnetic helicity without significantly perturbing the cross helicity. Runs R17a, R17b, and R17c use the above 
combined magnetic field, but also have a velocity that combines the Taylor-Green and ABC flows in a similar fashion: 

(£isinfc„ia; cos k v \y cos k v \z + £2 [B cos(k v2 y) + C sva.(k V 2z)]) , 
(— £1 cos k v \x sinfcuiy cos k v \z + £2 [A ain(k V 2x) + C cos(k V 2z)]) , 
£2 [Acos(k v2 x) + Bsm(k v2 y)] ■ 

In the case of the three R17 runs, both the magnetic field and velocity initial conditions have £1 = 0.93 and £2 = 0.07. 
The magnetic field is such that both portions of the initial condition are at k = 3, but the velocity has the Taylor- 
Green portion of the flow at k — 2 and the ABC portion at k = 3 initially. This allows for a perturbation in both 
the magnetic helicity and cross helicity. R20 involves a velocity which is a combination of the modified Taylor-Green 
velocity and a Beltrami ABC flow, such that: 

v tg'+abc = v tg'+abc ^ iCQskviX sinkviy C os k vX z + £ 2 [B cos(k v2 y) + C s'm(k v2 z)}) , 
v tg'+abc _ v t g '+ abc (~£ lS - m k vlX cosk vl y cos k vl z + £ 2 {Asm(k v2 x) + C cos(k v2 z)}) , 

V TG'+ABC = v TG'+ABC (2[Acos{kv2x)+BsHkv2y)] 

In run R20, £1 = 0.97, £2 = 0.03, and k v % and k v2 are set such that both the Taylor-Green and ABC porions of 
the flow are at k = 2. When combined with an ABC magnetic field this results in the addition of cross-correlation 
between the velocity and magnetic fields of 0.11 relative to the total energy. 

The type "£1 OT+£ 2 ABC" is an initial condition for the velocity which is the OT vortex and combines the OT 
and ABC magnetic fields, such that: 

b OT+ABC = &o^c ( ^ [_ 2asill 2fc ml y + asmfc ml z] + 6 [Bcos(fc ro2 y) + Csin(fc m2 z)]) , 
b OT+ABC = b OT+ABC [2 a sin fc m ix + asinfc m iz]+ £ 2 [Asm(k m2 x) + C cos(k m2 z)]) , 
b OT+ABc = feOT+^c^fasin^^ + asinfc^j + ^^cos^^+Bsin^ay)]) , 

where £1 and £2 set the relative fractions of OT and ABC respectively. Each of the flows of this type are such that 
k m i = k m2 = 1. 

Some runs were performed for more than 1000 tnl, where tnl = Lq/Uq is the turn-over time, and the maximum 
number of modes in the largest runs on grids of 256 3 points is in excess of one million. 



E. Global properties for all the runs 



We show in Fig. [I] the temporal evolution of the Reynolds number (left) and of the ratio of magnetic to kinetic 
energy (right) for most of the runs. The color table and symbols for runs is also given in Fig. [lj Since the runs 
are performed at relatively modest Reynolds numbers and numerical resolutions, but for long times, the Reynolds 
numbers eventually enter a regime of exponential decay where nonlincarities are weak. The burst of energy for run 
R20 at t ~ 250 is associated with the end of a plateau in the ratio Em / By and with a weak Lamb vector C = v x to 
(see Fig. |4]below). Examining the energy ratio, it is clear that two main regimes develop in these runs: some are close 
to equipartition, with a tendency to have an excess in magnetic energy as predicted by the statistical ensembles, and 
one where the magnetic energy wins all, and presumably under the influence of a strong relative magnetic helicity and 
an accumulation of Hm at the gravest mode of the computation. The run done on a grid of 32 3 points has its Reynolds 
number getting too low and as a result behaves considerably differently from the other Taylor-Green "A" flow runs 
on grids of 64 3 and 128 3 points. Of course v = is a possible solution of the MHD equations; this corresponds to the 
hydrodynamic attractor which can also be fluid in the forced case when the magnetic Reynolds number is too low. 
Also note that run R15, which is the Taylor-Green A flow perturbed by 1% ABC magnetic field, moves towards a 
kinetically dominated state; by the end of the run, this flow has kinetic energy dominating over magnetic energy by 
approximately a factor of 10 in the gravest mode. 



II. THE ROLE OF ACCURACY AND SYMMETRIES 

The ensemble of runs analyzed in this paper is shown in Fig. [2] in a plane first introduced in [HUE]; it delineates, 
in terms of the total energy, the relative importance of the two helical invariants. A peculiar feature of the I and A 
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FIG. 1: Left: Evolution of the Reynolds number denned as Re = VrmsLint/v where Lint is the integral scale (see text) as a 
function of time for most of the flows in Table [I] The plot is in lin-log coordinates and time is measured in units of turn-over 
time for each configuration. The inset shows a blow-up of the results for early times. Each run is indicated by symbols (and 
colors online) that are indicated in the center, a choice that is also followed in Figs. [2] and [4] . Right: Magnetic to kinetic 
energy ratio for most of the runs in Table |T] in lin-log coordinates. Many runs gather around quasi-equipartition, either above 
or below but for some runs, mostly those in which the relative amount of magnetic helicity grows significantly (such as runs 
R16 - R21b), this ratio gets substantially larger than unity. 



Taylor-Green runs is that, unless perturbed, they stay where they started, even though in these runs the symmetries 
are not imposed at all times. This may be related to the fact that it can be shown that, in the context of the fluid 
equations, symmetries are preserved by the dynamical evolution, a result that one may be able to extend to the MHD 
case [40]. Unperturbed, these two flows do not evolve in parameter space at these low Reynolds numbers. In the 
presence of perturbations, they do cover parameter space and evolve towards configurations with either strong Hm 
(and thus high ratio Em/Ev, as in the case of runs R16, R17a, and R17b), or strong He with near equipartition 
of kinetic and magnetic energy. Note that, in single precision and for long times, the I flow is perturbed by the 
accumulated round-off errors and it evolves toward another attractor, as shown in Fig. [3| whereas the accurate 
computation which maintains all symmetries evolves toward presumably a magnetically-dominated Taylor state, the 
errors introduced by insufficient precision lead to a quasi-equipartition of energy. This same behavior, where the 
single precision computation evolves towards a different attractor after sufficiently long times, is not evident in either 
the A or C Taylor-Green flows; however similar effects are observed when random noise is explicitly added to these 
initial conditions. 

The C flow also exhibits unique behavior as compared to the I and A flows in that, over the course of the computation, 
the accurate double precision run does not remain at the origin of the plane and instead moves to a state with 
2Hc/Et — —1 and Hm/Et = 0. As a result, the C flow without perturbations reaches an equipartitioned state, as 
opposed to a magnetically dominated state. 

We have also performed a more controlled and specific perturbation of the Taylor-Green symmetries by adding a 
fraction of a Beltrami ABC flow to the magnetic field and/or velocity of the Taylor-Green initial condition in the A 
configuration, i.e. by perturbing the flow explicitly with non-zero helicity (Runs 15-17c and 18-21b). By varying 
the amount of helical-ABC relative to non-helical Taylor-Green in the magnetic field, the value of Hm / Et can be 
adjusted in a controlled manner at t = 0. By also adding a fraction of ABC to the velocity, a set amount of 2Hq/Et 
can additionally be introduced to the flow. With larger perturbations to the magnetic helicity and cross helicity, such 
as in runs R16 - R17c, the symmetries are clearly broken and the runs reach the boundaries of the parameter space 
as predicted by the minimum energy principle (see Fig. |2]). Note that run R17c, which is performed on a grid of 
256 3 points, is only run for IOOtjvl. If this run were continued to longer times, as the other two R17 runs are, this 
run would likely reach the boundary. In the case of run R16, only the magnetic field initial condition is perturbed 
with 10% ABC, and the flow achieves a magnetically dominated state with nearly maximal Hm/Et- The three R17 
runs, which all have 7% ABC in both the velocity and magnetic field initial conditions, but are performed at different 
Reynolds numbers, evolve to a state on the boundary with both nonzero Hm/Et and 2Hc/Et when given enough 
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2H C /E T 

FIG. 2: Map of the trajectories taken through the (2Hc /Et , Hm/Et) parameter space by the various runs in Table [I] The 
color scheme is the same as that given in Fig. [I] with "x" marking the start position of each run and "o" marking the position 
at the final time of the run. Three runs (R5a, R5b, and R6) follow the black dotted line which goes from the origin to the point 
(-1,0) and 8 runs (Rla, Rib, RlOa, RlOb, RIOc, RIOd, Rlla, and Rllb) remain at the origin of the parameter space. Note 
that the origin of the parameter space is unstable and when symmetries are broken by perturbations the resulting trajectories 
can go in any direction depending of the phases of the perturbation. The energy minimization principle predicts that flows 
will move towards the boundaries of this space; that is 2Hc / Et = ±1 when the end value of Hm/Et < 0.5, and on the 
marked ellipsoidal curve bounding the upper portion of the plot when the end value of Hm/Et > 0.5. The parameter space is 
symmetric for negative values of Hm/Et- 



time. Note that although the three runs have the same initial conditions, they have take different paths through the 
parameter space with differing Reynolds number. 



III. THE INTERPLAY BETWEEN HELICAL INVARIANTS 



A. Magnetic helicity relative growth 



According to the equations written in \ I C the helical invariants play a central role in the evolution of MHD 
turbulence. Since the Taylor-Green flows have no helicity, we now examine a set of evolutions for several helical 
configurations that have been studied in the literature, namely the ABC (Beltrami) flows, the Orszag-Tang vortex 
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FIG. 3: Comparison for the magnetic to kinetic energy ratio for the Taylor-Green I flow between two codes, one implementing 
symmetries (run Rib, black circles), the other one not, (run Rla, black solid line), both in double precision, together with 
another run using single precision (run R2, blue dashed line). The double-precision in GHOST is sufficient to maintain the 
symmetries up to T=500, in units of turn-over times. Note the excellent agreement between the three computations performed 
on a grid of 64 3 points up to T ~ IOOtjvl, time after which the single-precision run departs from the others and evolves towards 
quasi-equipartition (at the final time, Em/Ev ~ 2, as often observed in the Solar Wind 15 ). The red dashed-dotted line and 
the green dotted line are the Taylor-Green I flow with added random noise of relative amplitude 10 -6 (run R3) and 10 -3 (run 
R4) respectively. Note that added noise behaves similarly to the single precision computation, but departs from the other runs 
and goes towards quasi-equipartition at earlier times. 



and some perturbations of such flows (see [ID and Table [T] for definitions). 

The Orszag-Tang vortex, without magnetic helicity, becomes highly correlated, but with the inclusion of some 
magnetic helicity, it evolves toward states which, as Hm/Et increases, are more and more magnetically dominated. 
With very small additions of magnetic helicity (runs R23a and R23b), Hm/Et grows to modest values at which 
2Hc/Et can still obtain a value of one and the flow has equipartition between kinetic and magnetic energy. However, 
with even a slightly larger addition of Hm (see runs R24 and R25 of Table the growth of H M /E T be gins to 
dominate and the flow moves towards more magnetically dominated states. Similarly, unperturbed ABC flows, with 
strong Hm, remain uncorrelated if initially so; but when perturbing them by adding some correlation between the 
velocity and the magnetic field, as in run R20, they follow similar evolutions but stay away from the singularity that 
occurs at maximum Hm/Et, He = 0. 



B. Vector alignment in MHD turbulence 



Relative alignment of dynamical fields (see eq. ^ for definitions) are shown in Fig. [4] for many runs, using the 
same line (color) encoding as in the preceding figures. A lot of runs reach an Alfvenic states (pc ~ ±1), some more 
slowly, and a few stay at low values: it has been known for a long time that the correlation between the velocity and 
the magnetic field grows with time (see e.g. |32j). 

Magnetic helicity seems more discriminating insofar as the long-time behavior of the runs: either k m i n HM / Et 
remains rather low, or else it approaches its maximal value. The runs that approach near maximal values of Hm/Et 
are those in which a sufficient amount of magnetic helicity is present in the initial condition. In Fig. 2, as well as in 
Fig. 4 it, can be seen that all the runs which attain large values of Hm/Et have at least slightly larger initial values 
of magnetic helicity than those that remain near zero. On the other hand, the normalized kinetic helicity shows a 
more varied set of behavior, with sometimes strong fluctuations between aligned and perpendicular fields, as for Run 
R20 (thick solid, red line) which evolves towards strong magnetic helicity (see the map in Fig. [2]). Note that Run R25 
(thick dash, purple, line) is likely to evolve in a similar manner. This analysis suggests that one ought to look more 
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FIG. 4: Temporal evolution of helicity in relative terms for all the runs of Table |l| with H C /E T (left), k min H M /E T (middle) 
and pv (right, see definition in eq. |3]l). The symbols and color scheme are the same as that in Fig. [I] Note the clear distinction 
between the runs for the two invariant quantities (either staying at values close to zero or close to their extrema), and the more 
varied evolution for pv (right), although its evolution is also influenced by magnetic helicity (see eq. Q). 
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FIG. 5: Probability distribution functions at T = IOtnl of the cosine of the angle between the magnetic potential and the 
magnetic field (top row) and the velocity and magnetic field (bottom row), with all PDFs normalized to their maximum value. 
Left-most: I flow runs Rla (black solid line), R2 (blue dashed line), and R3 (red dotted line); Middle left: C flow runs R5a 
(black solid line), R7 (blue dashed line), R8 (red dotted line), and R9a (green dashed-dotted line); Middle right: A flow runs 
RfOa (black solid line), R12 (blue dashed line), R14 (red dotted line), and Rf6 (green dashed-dotted line); Right-most: OT 
flow runs R22 (black solid line), R23a (blue dashed line), R24 (red dotted line), and R25 (green dashed-dotted line); see Table 
Ufor the nomenclature of the runs. 



in detail at the alignment properties of the various fields, by examining PdFs of the angle between various vectors. 
We show in Fig. [5] alignment probability distribution functions for several variables for several flows (see caption) , 
after 10 turn-over times: the top row is for magnetic potential and magnetic field, i.e. we are examining magnetic 
helicity, whereas the bottom row is for the velocity and magnetic field, i.e. we are concentrating then on the cross 
helicity. At T = 0, all undisturbed flows have a strong central peak corresponding to orthogonality of vectors (either 
a 1 b or v 1 b, and thus strong non-linearities, except for the OT case (right) for which pm peaks symmetrically at 
values slightly greater and slightly less than zero and pc indicates that there is a significant fraction of highly aligned 
velocity and magnetic field vectors. 

These PDFs confirm the results illustrated in Fig. [2] in showing an evolution towards either alignment of the 
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velocity and the magnetic field, or of the magnetic field and the potential, once the flows are perturbed, the more 
so the larger the perturbation. The I flows (left-most column) are in fact the harder to perturb insofar as alignment 
does not really develop and one only observes a widening of the PdFs around zero, i.e. a distribution of angles 
that remain nevertheless close to tt/2. For the C flow family of runs (middle left column), a perturbation at the 
level of 10~ x simply widens the distribution of angles (blue dashed curve), but increasing this perturbation leads to 
a totally different behavior and a flat distribution for magnetic helicity, the perturbations being less significant for 
cross-helicity in the sense that the PdF is changed but the overall distribution (its shape) is similar in all cases.. The 
A flows (middle-right column), as the flow is more perturbed from its highly symmetric initial conditions, the fields 
become more aligned, with an almost equal distribution for v and b (pc ~ ±1), whereas a clear alignment develops 
for a and b (pm ~ ±1). Finally, the Orszag-Tang flow (right-most column) starts from a different configuration of 
vectors, and its evolution as it is more perturbed is not so dramatically different (except for the green dash-dotted 
line which has a 50% OT-50% ABC mixture). 

On the other hand, the relative kinetic helicity (corresponding to alignment of velocity and vorticity) does not seem 
to follow a clear organization, unless magnetic helicity is strong and kinetic helicity does follow Hm and grows in 
relative terms, under the influence presumably of Alfven waves due to the large-scale magnetic field, as predicted in 
[37j . Growth of v— u or v— b alignment is a dynamical property of the Navier-Stokes or MHD equations, corresponding 
to the mutual interactions of shear and vorticity or shear and magnetic field [38]; in fact, such alignments properties 
have been found between all relevant fields, to different degrees [39] . This alignment property is shared by the rather 
low Reynolds number computations performed here, and for long times as well. 

C. Is there a dynamically significant ratio in MHD turbulence? 

We finally examine the relative role of the velocity and the magnetic field, in terms of energy distribution and 
timescales, for the runs performed at the highest Reynolds number, and thus with the largest extent of the inertial 
range (and numerical grid resolution). In order to do so, we look at the behavior in the inertial range around the 
peak of dissipation, when the turbulence is developed and the Reynolds number has not decreased substantially yet, 
of the ratio of magnetic to kinetic energy rs(fc) = EM(k)/ Ey(k) and of the eddy turn-over time to the Alfven time 
r T {k) = TNh{k) /TA(k), the latter being built on the magnetic field in the gravest mode. 

The ratio of timescales behaves as expected when evaluating the turn-over time on the energy spectrum variation 
with wavenumber. This can be seen in Fig. [6] for all flows, r T (k) increases with decreasing scale because of the way 
these two characteristic times change with scale, i.e. 1/fc for ta and [k 3 'Ev{k)\~ 1 / 2 for the eddy turnover time. On 
the other hand, and again for all flows, the energy ratio is constant and of order unity (but systematically slightly 
above in fact, as also regularly observed often in the Solar Wind), except in the largest scale in which it is dominated 
by initial conditions and the reinforcement of magnetic energy in the largest scale in all the runs dominated by an 
evolution towards the top of the map displayed in Fig. |2j and corresponding to cases with strong magnetic helicity. 
We also note that at later times of the order of three times the peak of dissipation, this result still holds but with, in 
all cases displayed here, an increase in Em/Ey at the largest excited scales, by a factor of 2 to 40 (not shown). All 
other runs of this study behaved similarly, as long as the Reynolds number is sufficiently high for turbulent mixing 
to take place. 

IV. CONCLUSIONS 

We have shown in this paper that the Taylor-Green configurations studied in |16) for their energetic properties, 
depart from their strong symmetries given a strong-enough perturbation, as can be encountered in high Reynolds 
number flows. They then evolve toward different characteristic behaviors in the Hm, Hq plane (map of Fig. |2j); these 
can depend on the perturbation (whether it has cross-correlation or magnetic helicity, leading the flows to different 
end-states in particular). It is not clear however if, when studying these flows at substantially higher Reynolds 
numbers, as was done in |16| but, contrary to [16j . not imposing the symmetries at all time, one will still have three 
different scaling laws for the total energy spectra for these three configurations. We also confirm that statistical 
mechanics, with an energy minimization principle, is an excellent predictor for the behavior of turbulent flows, as 
argued in p], and studied at moderate resolutions in [28], but the question still remains as to whether this leads to 
different scaling laws in MHD turbulence at high Reynolds numbers. 

Quasi-equipartition between the kinetic and magnetic energy is expected, on the basis of mixing of complex systems 
with a large number of degrees of freedom, although as shown in |29j , magnetic helicity, alone or in the presence of 
cross-helicity, may well prevent this from happening. What we have shown in this paper is that, in some cases 
with strong phase relationships such that the nonlinear terms are weakened considerably through alignment of the 
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FIG. 6: Nonlinear to Alfven time ratio (left) and magnetic to kinetic ratio (right), as a function of wavenumber at peak of 
dissipation for the following computations performed on grids of 256 3 points (see Table [i] for the nomenclature): run R9b (C 
flow with strong perturbation, black solid line), run 17c (A flow with 7% ABC, blue dashed line), run R21b (TG velocity and 
ABC magnetic field, red dashed-dotted line), and run R23b (99% OT and 1% ABC, green dotted line). Note the constancy 
of quasi-equipartition of energy throughout the inertial range, to be contrasted with the increase in the ratio of characteristic 
time scales in that same range. 



relevant fields (vorticity, velocity, magnetic field, magnetic current), other solutions are reachable with quite different 
properties. This is a bit akin to a potential flow in hydrodynamics: when the vorticity is identically zero, this is an 
exact solution of the Navier-Stokes equations although an unstable one, and vorticity, like a seed magnetic field, grows 
over time. The reason why these solutions do not destabilize in a Lyapounov time, which can be close to an eddy 
turn-over time, is probably due to the fact that symmetries are very strong properties of flows that are preserved by 
the dynamical evolution |40[ 141 J ; would stronger turbulent flows at higher Reynolds number lead to more complex 
behavior than what is observed here? Or, since MHD turbulence is known to be a dissipative system including in the 
limit of infinite Reynolds number , is the emergence of large-scale helical coherent structures enough to insure a 
simple evolution following an energy minimization principle? 
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